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Abstract 

The model with the meson spectrum, consisting of zero-width equidistant resonances, is 
considered with connection to current correlators in coordinate space. The comparison of the 
explicit expressions for the correlators, obtained in this model, with the experimental data of 
ALEPH collaboration on r-decay is made and good agreement is found. 



1 Introduction 



The main objects of our studies are correlators of charged hadron currents with the value of isospin 
/ = 1 

p!fJi^) = {jl'\x)ji^\oy) (1) 

for vector and axial channels, i. e. for 

j;^\x) = uj,d; jI^\x) = uj,j,d. (2) 

These objects play important role in QCD, since they contain information of two kinds: one 
is the spectral features of hadrons - the set of the mass poles and the corresponding quark 
coupling constants c„, which define behavior of the correlators for large x, on the other hand, at 
small distances, where perturbation theory works, and OPE is assumed to be valid, the coefficients 
of the correlation functions expansion give us matrix elements of local field operators - the vacuum 
condensates like (G^,^). 

Current correlators were studied long ago in the framework of perturbation theory, which 
works at small x because of asymptotic freedom of QCD [l]. Power corrections to this behavior, 
expressed in terms of gluon and quark condensates, were found with the help of OPE, which is 
consistent also only for small enough distances [2]. There are many papers, where calculations 
of ^ are given on the base of different approaches. So, the random instanton liquid model, 
where it is supposed, that the nonperturbative QCD-vacuum is dominated by strong classical field 
configuration - instantons, whose position and color orientations are distributed randomly, - gives 
good result for the difference of the correlation functions in all region of distances (0 < x < 1.5 fm) 
[H]. Correlation functions and their moments in momentum space are considered in some papers 
with the purpose of extraction from them the values of the lowest OPE condensates f5j, [H]. 
There is another model - the so-called Minimal Hadronic Ansatz (MHA), in which the spectral 
density consists of the pion state, a vector state and an axial state plus continuum. In this this 
approximation one can compute corresponding condensates [I| or compare the moments of the 
polarization operators with the experimental data |H]. The method of calculation of correlators, 
more general, than MHA, was proposed in the paper 0. This method, as well as MHA, assumes 
the limit of large number of colors Nc, when the hadronic spectrum is represented by the set of 
zero- width levels. One can calculate the spectrum quasiclassically [in], restore the correlator in 
momentum space by this spectrum and by residues and then go to coordinate space. In the same 
paper was made the comparison of obtained expressions for the correlators with the lattice data 
[TT] and [12] and good agreement of the results was found. 

This paper continues the line of research, started in jHj. We will be interested in explicit 
expressions for current correlators in the model with equidistant spectrum, and the main 
purpose is to compare these expressions with the real experimental data of ALEPH collaboration 
[13] on hadronic r-decay. Note also, that the model of zero-width equidistant resonances has been 
considered with respect to correlation function in [H], jTH], where the main attention was spared 
to the momentum-space representation and the corresponding OPE coefficients were estimated. 
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We will be interested in the current correlators, calculated in the Euclidean coordinate space and 
compare them directly with the experimental curves. 



2 Current correlator: the model and experiment 

We work with QCD in the limit of large number of colors A''^. In the framework of this limit 
the hadronic spectrum is a set of zero-width resonances. Masses m„ and residues c„ of current 
correlators are stable at large Nc limit and define the imaginary part of the polarization operators 
(see Appendix): 

^ oo 

-lmll{s) = CnS{s - ml), (3) 

TT ^ — ^ 

n=0 

On the other hand, one can obtain the equidistant mass spectrum and the set of residues in 
analytic form for all n quasiclassically, solving the corresponding Hamiltonian problem for two 
quarks 01101: 

where ml are the masses of low lying resonances in each channel, and the quantity = Area is 
defined universally for all channels, where a is physical string tension. 

In that way one can easily find with the help of the correlation function in the momentum 
space and then go to the coordinate one (see Appendix). 

It is convenient to consider the normalized correlation functions and to define the ratio R^'^^x) 

as 

R^^\x) = ^^^^^^^ ^ 9ZftM^ (5) 

where 

6 

Pfreeix) = (6) 

Experimental data provide us with information about the spectral functions Vi{q'^) and ai(g^) 
in the interval < s < Sq = 3.5 GeV^: 

Imn('^)(s) = ^v^{s), lmU,{s) = ^a,{s), (7) 

i. e. about the imaginary parts of polarization operators Il^'"\q'^) and ni(g^), defining the correla- 
tion functions in momentum space (fT6jl and (fT7|l . 
Using the dispersion relation 

oo 

.2^ /• 1 1 



n(g^) = / dt- --ImUit) (8) 

t — q"^ TT 
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and the known representation for the McDonald function Ki{z) for z = y/—sx'^, it is easy to obtain 
the spectral representation for correlation functions in coordinate space: 

oo 

P^'\xe) = J sp^'\s)D{yG,XE)ds, (9) 
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where xe = V—x^ is a distance in Euclidean space, p^'"'^"'\s) = -z—:^vi{s) \ai(s) are corresponding 
spectral functions, and 

D(v^, xe) = -^^Ki{^sxe). (10) 

One should take care of the term q^qyli2{,<f') from (fTTjl . which corresponds to the contribution 
of the pion state. It is taken into account by adding to P^"-\xe) the term (see, e. g. [^) 

Pi''\xE) = -f^mlD{m^,XE). (11) 
The integrals in considered combinations 

oo 

Pv±a{xE) = j s p^±a{s)D{^,XE)dsT flmlD{m^,XE), (12) 



where pv±a{s) = p^'"\s) ± p^"'\s), Pv±a = -P*-^'' '^P'^°'\ were taken numerically by Simpson method, 
using the experimental data for difference and sum of the spectral functions in the region < s < Sq. 
The spectral functions Pv±a{,s) were replaced by the parton model prediction for s > sq, i. e. 
fi(s) + ai(s) = 1, vi{s) — ai{s) = for s > sq. Under this replacement we do not take into 
account perturbation theory corrections, which change little behavior of the sum of the current 
correlators for small x. However, this perturbative corrections, corresponding to the hybrid states 
in the spectral sum in terms of our model, are not taken into account also, so the disregarding of 
them in the experimental curves is coordinated with the model calculations. 

The experimental curves, normalized to the free answer and the corresponding model 
predictions are shown on the Fig. [T] for two values of the cut-off scale Sq = 2.3 GeV^ and 
So = 3.1 GeV^. Notice, that the choice of the cut-off parameter Sq changes noticeably behav- 
ior of Pyj^ai^E) only in the intermediate region of x, pushing away the theoretical curve from the 
experimental one. In addition, the error bars increase with increasing sq in this region because of 
large uncertainties of the experimental data for sq > 2.8 GeV^, so agreement of the experimental 
curves with the model predictions remains acceptable even on the quantitative level. 

On the other hand, using meson masses and decay constants as fitting parameters, one can 
place the curves, obtained in the framework of the model at large Nc, inside the error corridors 
completely (see Fig. [2l). The corresponding values 

ml = 0.62 GeV^ Aj = 0.048 GeV^ 

= 135 MeV, U = 89 MeV, (13) 

m„„ =1.18 GeV 
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reproduce the table values (21) with 2 — 4% accuracy. 

The errors was computed with the help of experimental covariance matrices 

Pv±a{s., S') = (Ap„±„(s) Apy±a{s')) (14) 

by the standard expression 



3 Conclusion 

We have considered the model, based on QCD in the limit of large number of colors N^, in 
which the meson spectrum is taken as the series of zero-width equidistant levels. This model is 
surely far from the real physical world and is over simplified. However, as can be seen from the 
results of comparison of the calculations, performed in the framework of this model, with the high 
precision experimental data of ALEPH collaboration for the correlators of hadron currents, the 
knowledge of several spectral features (in our case, these are the data about low lying resonances 
in corresponding channels and the correct behavior of the mass spectrum and of the set of residues 
for n — > oo) is enough to reproduce correctly on qualitative and even on quantitative level such 
integral objects (in a sense of the spectral density) as correlation functions (P). Besides, note also, 
that all the parameters of the model (meson masses, decay constants, string tension) are fixed 
by their physical values, and calculated correlation function have in principle a form of a fixed 
prediction without any fitting parameters, however, having shifted quite a little the values of this 
input parameters (~ 2 — 3%), one can reach a very good agreement with the experimental data. 
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so so 

y (is y ds' Pv±a{s-, s') S D{^, Xe) s D{Vs', Xe) 



1/Z 



(15) 
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Appendix 

For conserved vector current one can write the following expression in momentum space: 

P^':){q) = ^ I d'xPi:){x) exp(zgx) = {q,q, - g ,,)Il^'"\q^) . (16) 



There are two Lorentz structure in the axial channel: 



Pl:){q) = ^ j d'xP^:\x) expiiqx) = {q^q, - q''g^,)Ur{q^) + q^q^H^iq'). (17) 

The normalized correlation functions in the Euclidean coordinate space are given by following 
expressions in accordance with result of ^: 

R^^\z,) = izlK,{z,) + ^-^p{K, z,) (18) 

R^''\^a) = -^/V^f^l(7^a) + ^P(foa,^a), (19) 



which are depend on the dimensionless coordinates Zc = v— a?mc, where rric are the masses of the 
resonances in vector (p - meson) and axial (oi - meson) channels. The contributions of the pion 
and the p - meson states are taken into account separately and are described by the first members 
in R^^\x). The parameters 

o \mfj 477^/ nia rUa 

are fixed by the physical values of meson masses, decay constants and string tension constant: 

ml = 0.6 GeV^ Aj = 0.047 GeV^ 

= 140 MeV, U = 92 MeV, (21) 
rUa, = 1.23 GeV, a = 0.17 GeV^. 

Remaining terms, containing the function 

oo 

v(b z) - f- exp (-- - u] 1 + ^^P(-H(fe-1) .22) 
J u'^'^^y Au (l-exp(-H)' ' ^ ^ 



are determined by contribution of a set of the equidistant levels with the masses 

{m^^^Y = rn^n + ml, m^ = Ana. (23) 
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Fig. 1. Normalized Euclidean coordinate space correlation functions Pv±a{x) for sq = 2.3 GeV^ 
and So = 3.1 GeV^ 
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Fig. 2. Normalized Euclidean coordinate space correlation functions Pv±a{^) for the set of param- 
eters (13) and Sq = 2.5 GeV^ 
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